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ABSTRACT 


In this thesis, we have studied two different problems in the area of stability of fluid flow 
past a deformable solid. First we studied the stability of Newtonian fluid flow past a 
deformable wall using a combination of asymptotic and numerical methods for the case 
of Couette flow as well as Poiseuille flow past a deformable wall. The fluid has viscosity 
p thickness R supported by a viscoelastic medium of elasticity G and viscosity pg and 
having thickness H R .We have shown that 'wall mode instability' predicted by the 
asymptotic analysis is also captured by numerical methods for both Couette flow as well 
as Poiseuille flow past a deformable wall. Although in Couette flow the coupling 
between mean flow and fluctuations is only via an additional term in the boundary 
condition for the tangential velocity at the interface, as compared to Poiseuille flow 
where other than this one extra term in Orr-Sommerfeld equation which is responsible for 
the coupling. Further we have shown that there exists a mode in Couette flow which has 
scaling as F ~0 (Re'^) ,(r = p if/GR, Re - U*R p/ \i) as compare to wall modes where F 
~0 (Re'^ and further we confirm that this mode is also present in the case of Poiseuille 
flow. To examine the nature of the modes, we calculated the scaling of boundary layer 
thickness, velocity and velocity derivatives as a function of Reynolds number (Re) for 
both the cases and we find that new modes are indeed inviscid modes. Further we find the 
boundary layer thickness as a function of Reynolds number and this also indicates that 
the new mode is an inviscid mode. Other than these, the distribution of the perturbation 
velocity field also confirms validity of our results. 

The second problem we studied concerns the stability of two-layer Newtonian 
fluids flow past a deformable for the case of plane Couette flow. The fluids have 
viscosities pi, P2, thickness P R, (1 - P) and densities pi, p 2 respectively supported by a 
viscoelastic medium of elasticity G and viscosity p^and having thickness HR. Here, the 
idea is to explore the possibility of using deformable solid layers to suppress two-fluid 
interfacial instability due to viscosity stratification. We have shown that deformable wall 
could have both stabilizing and destabilizing effect on the two-layer flow past a 
deformable solid wall, depending on the arrangement of the two fluids and their relative 
thicknesses. Further we have shown that deformable wall has an effect on long and 



intermediate waves but it does not have any effect on short wave perturbations. Since 
short wave perturbations are confined near the interfaces, therefore there is no possibility 
of interference in the case of short waves. However, by ensuring sufficient interfacial 
tension between the two fluids, it is shown that it is possible to suppress two-fluid 
interfacial instabilities at all wavelengths by the deformable solid layer. 
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CHAPTER 1 

1.1 Introduction 

Fluid flows are ubiquitous in the technological and the natural world. Their importance 
and abundant use in science and technology result in a long list of areas that are used 
every day. Recently, one of the important areas of study in fluid mechanics is fluid flow 
past a deformable solid wall. The stability of fluid flow near a deformable wall is of 
interest in diverse applications such as marine and aerospace propulsion, biotechnology 
and polymer processing applications. The flow of blood and other fluids in the body takes 
place through deformable tubes, and the separation and purification processes in bio- 
technological processes often involve flow in a tube and channels made up of polymeric 
materials. Apart than these situations, it is also very important in polymer processing 
applications such as multilayer co-extrusion, where there is need to control flow in such a 
way that there is no possibility of mixing between the layers, so as to enhance the product 
quality. 

From a technological viewpoint, the presence of instability has two effects: it 
increases the drag force and the power consumption, and it also enhances mixing and 
heat and mass transfer. In processes where it is advantageous to have lower drag force, 
the wall and fluid properties could be adjusted so that the flow is stable and remains in 
the laminar regime. Alternatively, in situations where greater mixing and heat or mass 
transfer are desirable, the fluid and wall properties could be suitably chosen to make the 
flow unstable. Therefore, a detailed imderstanding of the mechanism(s) responsible for 
destabilizing the flow could be used in biotechnology applications. Other than this, we 
know fi'om the past studies of Yih [1967] and Hooper & Boyd [1983] that multilayer 
flows of Newtonian fluids are unstable even in the limit of small (but finite) Reynolds 
number. Thus by delaying the instability of multi-layer fluid flows at finite Reynolds 
number, we can reduce the mixing during multi-layer extrusion of polymers. For that 
purpose, we explore the possibility of using a deformable wall and we will show further 
in Chapter 3 that it has the ability to suppress the interfacial instability. Thus, an 
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understanding of the stability of flow past a deformable wall is often crucial in order to 
obtain the desired results. 

There is a significant difference between instability criteria for the fluid flow in a 
rigid channel and fluid flow through deformable wall. In a rigid channel for the high 
Reynolds number flow, the fluid inertia plays an important role in destabilizing the flow. 
However, the inertial effects are absent in a low Reynolds number flow, and instabilities 
have to be driven by a different mechanism in this case. The low Re 3 molds number flow 
of a Newtonian fluid past a rigid surface is always stable because the Stokes equations are 
not explicitly dependent on time, but the flow of viscoelastic liquids could become 
unstable due to the time dependence of the elastic terms in the momentum equations. 
When a fluid flows past a deformable wall, the dynamics of elastic wall and fluid get 
coupled and cause an additional time dependence via the deformable wall, which causes 
instability at very low Reynolds number as shown in the experiment of Krindel & 
Silberberg [1979]. Hence an anomalous drag force is observed when the fluid velocity is 
increased beyond a critical value, which depends on the fluid viscosity and the wall 
elasticity, in addition to the usual dimensionless groups characterizing the fluid domain. 
The nature of the instability observed is also different in a deformable wall channel; there 
is a gradual increase in the drag force as the Reynolds number is increased, in the contrast 
to the sharp (nearly discontinuous) jump in the drag force in the rigid channel at the 
critical Reynolds number. In the case of two-layer fluid flow in a rigid channel, the flow 
is unstable even in the range of very small Reynolds number but we will show in our 
results in Chapter 3 that it is possible to stabilize this type of flow by making solid layer 
deformable. 

The study of high Reynolds number flow of a fluid past a deformable wall has 
been motivated by the desirability of drag reduction in marine and aerospace propulsion. 
There have been many papers after the first significant study done by Kramer [1957, 
I960]. In this thesis, we consider the stability of single layer fluid past a deformable wall 
at modetarate and high Reynolds number, and this is an extension of the work of 
Kumaran, Shankar, Srivatsan and Carpenter for flow past a deformable solid wall (their 
works are discussed in an appropriate place). Here we briefly discuss their previous work 
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related to our present work. The stability of Couette flow past a deformable wall first 
studied by Kumaran et.al [1994] in the zero Reynolds number limit and they found that 
flow is unstable even for vanishing Reynolds number. Later this zero Reynolds number 
viscous mode was numerically continued to intermediate and high Reynolds number by 
Srivatsan and Kumaran [1997]. They found that this mode of instability is present for 
intermediate and high Reynolds number. Shankar and Kumaran [2002] studied Couette 
flow past a deformable wall, but they carried out an asymptotic analysis at high Reynolds 
number for wall modes. They refer the above modes as ‘wall modes’ since the vorticity 
for this mode is confined only in a region of thickness <5 near the fluid-solid interface, 
here <5 ~ 0{Re^^^) is the boundary layer thickness, which is small compared to thickness 
of the fluid. In that study they found that there are many modes of instability 
corresponding to the same class of wall modes. Further they continued the lowest mode 
(i.e. wall mode with the lowest numerical value of Cr) numerically to intermediate and 
low Reynolds number. While doing so, they found that the first wall mode is the same 
mode as that of Srivatsan and Kumaran who continued the low Reynolds number 
instability of Kumaran et.al [1994]. Before these studies, Kumaran [1995] predicted that 
inviscid mode is always stable for plane Couette flow past a deformable wall. Kumaran 
referred to the above modes as ‘inviscid mode’ because boundary layer thickness d ~ 
0(Re' ). Apart from the study of single layer Couette flow past a deformable wall, we 
also study the stability of Poiseuille flow past a deformable wall in a channel in order to 
examine whether all the modes found in Couette flow are present in Poiseuille flow or 
not. This problem had not studied extensively other than the study of Davis and Carpenter 
[1997], but they did not discuss all class of modes. Besides that motivation, the stability 
of Poiseuille flow is significantly different from that of Couette flow since it has 
singularities at high Reynolds number, which gives rise to a critical layer, similar to 
Shankar and Kumaran [1999, 2000]. 

The other problem we intend to study in the thesis is the stability of two layer 
flows past a deformable wall. The first significant study had been done by Yih [1967] for 
the stability analysis of two-layer flow. Yih did an asymptotic analysis of the interfacial 
mode for long waves and shows that jump in viscosity across the interface can lead to 
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long wave instability, even in the low Re limit. Later Hooper and Boyd [1983] solved the 
linear stability problem for two-dimensional plane Couette flow. They showed that this 
flow is unstable for short waves, and there is no instability in the long-wave limit. The 
asymptotic expression for the growth rate of the interfacial mode for short wave is 
essentially a sum of three terms. The term associated with surface tension, which is 
always stabilizing. The next term depends on gravity and the ratio of the densities and 
third term depends on ratio of viscosities: the last term is always destabilizing if the 
densities are equal. 

We now give an overview of the present thesis. In this thesis, we first intend to 
study the problem of single layer Couette as well as Poiseuille flow past a deformable 
wall in Chapter 2. For these problems, we would like to determine the nature of 
instability at intermediate and low Reynolds number corresponding to the wall modes 
present at high Reynolds number. Other than this, we are also interested to observe the 
effect of singularities present is Poiseuille flow past a deformable wall. 

For the above mentioned purpose, we first wrote the governing equation for mass 
and momentum conservation with the required boundary conditions. Then we non- 
dimensionalise the governing equations by using dimensional quantities: thickness of 
fluid, velocity of upper plate, viscosity of the fluid and density of deformable solid wall in 
such that the non-dimensional velocity gradient of fluid at fluid-deformable wall interface 
remains same. By doing so we ensure that boundary condition and Orr-Sommerfeld 
equations are equally valid for Couette flow as well as Poiseuille flow. For solving Orr- 
Sommerfeld equation numerically we used the principle of expansion of the solution as a 
linear combination of orthogonal solutions. But to solve Orr-Sommerfeld equation by 
mentioned method we need to know sufficiently close initial guess. Thus we perform 
asymptotic analysis for high Reynolds number limit and the results of asymptotic analysis 
we used as an initial guess in the numericeil code. First of all we calculated and presented 
the result of Shankar and Kumaran [2002] at high Reynolds number to verify the code. 
Further we calculated neutral stability curve and found a new mode. After examination of 
velocity field for the fluid, displacement field for deformable wall and scaling of all the 
quantities for the fluid as well as deformable wall. We concluded that this new mode is 
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nothing other than inviscid mode. Finally we tried for Poiseuille flow and found that all 
modes those are present in Couette flow is also in Poiseuille flow. 

Then further we studied the stability of two-layer Newtonian fluids past a 
deformable wall in Chapter 3. For this problem we axe first looking for constant density 
ratio equal to one and further we will show the effect on density difference on the 
stability. Hence out of the three terms (mentioned above) for short wave as shown by 
Hooper and Boyd only first and third are valid for equal fluid densities as mentioned by 
Joseph and Renardy [1993]. The case of equal densities of the fluid there exist some rules 
of thumb, which tells that we can stabilize long waves by placing the less viscous fluid in 
a thin layer, and stabilize the short waves by ensuring sufficiently high surface tension. 

In present study our aim is to show that: 

I. Two layer flows past a rigid wall is imstable and it is possible to stabilize with 
the help of deformable wall, 

II. Two layer flow past a rigid wall in stable and it is possible to destabilize the 
flow with help of deformable wall and finally, 

III. Effect of various parameters on the stability of our system of flows. 

For this analysis all the governing equation, non-dimensionalisation and numerical 
methods are same as in Chapter 2. But here we perform long wave and low Reynolds 
number asymptotic analysis as compare to high Reynolds number asymptotic analysis for 
single layer flow to find initial guess of numerical code. To check to validity of numerical 
code we recaptured the results for rigid wall case. Then first, we discuss the case where 
fluid flow in rigid channel is unstable (/i,> 1 & P > 1/2 or 1 & P < 1/2) and 
deformable wall has stabilizing effect on the flow. Again we discuss the case where fluid 
flow in a rigid channel is stable (//^> 1 & p < 1/2 or > 1 & p < 1/2) and deformable 

wall has destabilizing effect on the system of flow. Although from above cases we are 
able to know that deformable wall has stabilizing as well as destabilizing ability. But to 
calculate the range of elastic properties of deformable wall we need to know the neutral 
stability curve between gamma {T = /x U /G R) vs. k. Thus we plotted neutral stability 
curve for both the above cases. In this Chapter we finally discuss the effect of various 
quantities on the stability of two-layer fluid flow. 
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CHAPTER 2 


Linear stability analysis of high Reynolds number Poiseuille 
flow and Couette flow past a deformable wall. 

2.0 Introduction 

We first started with an asymptotic analysis in the limit of high Reynolds number to find 
the instability for Couette flow past a deformable wall, which is equally valid for 
Poiseuille flow. This is because, for both the cases non-dimensionalisation had been done 
in such a way that velocity gradients are same near the wall (V. Shankar, personal 
communication). The stability of Couette flow past a deformable wall first studied by 
Kumaran et.al (1994) for zero Reynolds number limit and they found that flow is 
unstable even for vanishing Reynolds number. Later this, zero Reynolds number viscous 
mode was numerically continued to intermediate and high Reynolds number by Srivatsan 
and Kumaran (1997). They found that this mode of instability is present for intermediate 
and high Reynolds number. Shankar and Kumaran (2002) studied Couette flow past a 
deformable wall, but this time they did asymptotic analysis at high Reynolds number for 
wall modes. They refer the above modes as wall modes since the vorticity for this mode 
is confined only in a region of thickness Here 5-0 (Re'^''^) is the boundary layer 
thickness, which is small compared to thickness of the fluid and is confined near the 
fluid-solid interface. In that study they found that there are many modes of instability 
corresponding to the same class of wall modes. Further they continued the lowest mode 
(i.e. wall mode with the lowest numerical value of Cr) numerically to intermediate and 
low Reynolds number. While doing so, they found that first wall mode is the same mode 
as that of Srivatsan and Kumaran when they continued the low Reynolds number 
instability of Kumaran et.al (1994). Before these studies, Kumaran (1995) predicted that 
inviscid mode is always stable for plane Couette flow past a deformable wall. They called 
above modes inviscid mode because boundary layer thickness <5 ~ 0 {Re ). 

In this present study, our first aim is to know the nature of neutral stability curves at 
intermediate as well as low Reynolds number corresponding to wall modes for Couette 
flow. Then further we intend to study the effect of base velocity profile on the Couette 
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flow stability modes by changing base velocity profile to Poiseuille flow. Further we 
would like to check critical layer thickness as suggested by Orr-Sommerfeld equation 
since second derivative in the case of Poiseuille base velocity profile is non zero, which 
reveals the possibility of critical layer singularity at high Reynolds number. Finally, we 
would like to explain nature and causes of the entire class of modes those are present in 
above plane Couette and Poiseuille flow past a deformable wall. Besides these things, it 
is important to study Poiseuille flow past a deformable wall in a channel on its own right. 
Since there is little information available in literature, mainly by Davis and Carpenter 
(1997) but they did not study the problem extensively. 

2.1 Problem formulation 

The problem consists of fluid layer of thickness R with base velocity profile as 
plane Couette and Poiseuille flow, base is supported by a viscoelastic solid of thickness H 
and upper layer is bounded by a rigid wall. The fluid is assumed to be an incompressible, 
Newtonian fluid with constant viscosity and interfacial tension. Here our aim is to 
examine the wall mode instability for plane Couette and Poiseuille channel flow and to 
determine the effect of base velocity profile on wall mode instability. Although we 
examine instability for both viscous and inviscid limit, we extend our solution to 
moderate Reynolds number using a numerical method with the asymptotic solution as 
initial guess. Governing equations for both fluid and viscoelastic solid for mass and 
momentum transfer in the bulk phases are the continuity equations and the Navier-Stokes 
equation. 



Rigid wall 


Fig 2.1 Schematic diagram of plane Couette flow past a deformable wall 
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Fig 2.2 Schematic diagram of plane Poiseiiiile flow past a deformable wall 

Z2 Governing equations 

2Z1 Governing equations in Dimensional form 

This system consists of a single layer of fluid flowing past deformable solid as shown in 
Fig 2.1 and 2.2. The governing mass and momentum conservation equations with 
required boundary conditions can be written as follows. 

Mass balance for fluid; 

a,v,. =0, 

Mass balance for deformable wall: 

=0 , 

The indices i, j = x and z. 

The momentum balances in the fluid and in the deformable wall are given by: 


( 2 . 1 ) 

( 2 . 2 ) 


p(.d,+ Vjd^)v^ = -d,p + /xd]v,. 

(2.3) 

/7(a>, + iG + M^d,)d]u, . 

(2.4) 
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F^ere we have assumed the density of both fluid and deformable wall to be the same for 
the sake of simplicity. Here repeated indices are summed over and single index means 
component of a vector. Here v, it, and pg represents velocity of fluid, deformation in 

viscoelastic wall, pressure inside fluid and pressure inside deformable wall respectively, 
p, p, and G represents density of fluid and deformable wall, viscosity of fluid, 

viscosity and coefficient of elasticity of the deformable wall respectively. 

"Now for solving equations from 2.1 to 2.4 we need boundary conditions at upper wall 
and deformable wall-solid interface and matching conditions at fluid-deformable wall 
interface. 

The no slip condition at rigid solid and fluid interface is (at z = R): 

V, =0, (2.5) 

The impermeability condition at the solid surface is (atz = R): 

V, = 0, (2.6) 


The tangential and normal velocity balance at fluid-deformable wall interface: 

v,=d,u^, (2.7) 


- 

The shear stress balance at fluid-deformable wall interface can be given as: 

The normal stress balance at fluid-deformable wall interface can be given as: 

/ 7 . Uj = n, a, , rtj , 


( 2 . 8 ) 


(2.9) 

( 2 . 10 ) 


The no slip condition at rigid solid surface is (at z = -HR): 

«, = 0 , ( 2 . 11 ) 

The impermeability condition at solid surface is (atz = - HR): 

( 2 . 12 ) 
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For fluid stress can be represented written as: 
Ty = -Pf fj. (d^Vj + djV.) , 


For deformable wall stress can be represented written as 
= -Pg ^iJ + . 

The dimensional stress ‘ Cy ’ consists term ‘ G (d^Uj + djU^) ’ the elastic stress due 
to the strain in an incompressible deformable wall and the terms ‘ p^d^id^Uj + djU^) ’ is 
the viscous stress due to the gradients in the velocity in deformable wall. 


2.2.2 Governing equation in non-dimensional form 


Before starting non-dimensionalisation we want to mention that in above section namely 
in 2.1 even though we did not put * anywhere, all the expressions represents dimensional 
quantities. But in this section dimensional quantities are represented by * superscript and 
non-dimensional quantities are represented without * superscript. The dimensional 
equations from 2.1 to 2.12 can be non-dimensionalised by using the following 
expressions: 

x=x*/R, z = z*/R, 

t = t* U^/R, 


V = V */U*, 

P = P*/(p U*/R), 
Re = U*Rp/p, 
Tij=T\/(plf/R).. 


Pr =Pg/p, 

T = p U*/GR, 

I = pGR^/p\ 
aij = a;/(pU*/R). 


Where x, z, t and v has obvious meanings as mentioned above. Here we represent 
dimensional velocity ’if' for non-dimensionalisation of other quantities. Where U*= Vmax 
for Couette flow and if = 4vmax for Poiseuille flow. We choose above method for non- 
dimensionalisation velocity in order to keep the non dimensional gradient of velocity at 
the interface of fluid and deformable wall to be the same for both flows. Here Vmca 
represents maximum velocity for the case of Couette flow as well as Poiseuille flow, F is 
the non-dimensional stain rate and Re represents Reynolds number. 
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Thus non-dimensional equations con-esponding to equations 2.1 to 2.12 can be 


written as follows. 

Mass balance for fluid: 

=0, (2.13) 

Mass balance for deformable wall: 

d,u.=Q , (2.14) 

Momentum balance for fluid and defonnable wall can be written as: 

Re {d, + V .a ^ )v,, = -f- d)v. . (2.1 5) 

Re(aV', ) = (2.16) 

The no slip condition at rigid solid and fluid interface is (at z =1): 

v,= 0, (2.17) 

The impermeability condition at the solid surface is (at z = 1): 

V. = 0, (2.18) 


The tangential and normal velocity balance at the fluid-deformable wall interface: 


V- = , 

The shear stress balance at fluid- deformable wall interface: 


The normal stress balance at fluid- deformable wall interface: 
n. T: n. = n. <7,,. tij , 

The no slip condition at rigid solid surface is (atz =-H ): 


(2.19) 

( 2 . 20 ) 

( 2 . 21 ) 

( 2 . 22 ) 

(2.23) 
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The impermeability condition at solid surface is (atz =-// ): 

=0 > (2-24) 

Non-dimensional stresses in liquid and deformable wall can be given respectively as 
follows: 

Tij=-Pf5,j+ 

^// = ^ Mrd,){djU,+d,u.). 

ZZSBase state calculation for both Poiseuille and Couette flow 

After solving two-dimensional Navier-Stokes equation in non-dimensional form for fig. 

2. 1 we are able to get base velocity profile as follows; 

V. = 

K =0, 

Base state deformation for the deformable solid in non-dimensional form is: 

u, = r(z-h//), 

= 0 . 

After solving two-dimensional Navier-Stoke equation in non-dimensional form for fig. 

2.2 we are able to get base velocity profile as follows: 

^ = z(^-z), 

T, =0, 

Base state deformation for the deformable solid in non-dimensional form is: 

= r (z + /f ), = 0. 
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23 Linear stability analyses 

23,1 General equation for linear stability analysis 



Rigid wall 

Fig 2.3 Schematic diagram of plane Couette flow past a deformable wall with 
perturbation. 


Linear .stability analysis predicts the length scale and the sufficient condition for 
instability. The objective of this analysis is to determine the stability of the fluid motion 
by applying infinitesimal perturbations to basic laminar flows. For simplicity we assume 
two-dimensional perturbations. 

The base states for Couette flow as well as Poiseuille flow are given in section 2.2.3. So 
now we add perturbations to the base profiles to find instability of flow. 

Thus perturbed quantities can be represented as follows; 

{v,-v^,v^-v^,u-u,Pf-p,,p^- p^] = (x-c/)}- 

Where, c and k are the growth velocity and the wave number of the instability 

respectively and base state condition is given in section 2.2.3. 

Applying this perturbation state to equation 2.13 to 2.24 and equating linear term on both 
the sides, we obtain: 

d,v^ + ikv^ = Q, (2-25) 
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(2.26) 


Re{ik (v,-c)v^+ «vjv,} = -ikP^+ 

Re{ik ( -c)vj = -d^P^+[d^^-k^] . (2.27) 

For cleforinable wall continuity equation, x and z component momentum balance for 


perturbed states are respectively; 


d.u^ + iku^=0. 

(2.28) 

- Re k^c^u^ = -ik + (^-?7,. ikc)[dl-k^]u^. 

(2:29) 

-Re k^c^u^ = -dj)^ + kc)[dl-k^]u^. 

(2.30). 

After combining above equations 2.25 to 2.27 for fluid and 2.28 to 2.30 for deformable 

wall we get a standard equation, which is known as Orr- Sommerfeld equation for fluid 

and deformable wall respectively. 

For deformable wall: 

-Re ec\dl-k^]u^= {~rijkc)[dl-efu.. 

(2.31) 

(l/('^)){ Re A:Vw^+ {~r]ri^c)[d'l-k^]u^ } . 

(2.32) 

For fluid: 


R^ik {(v, -c) [dl-k^]+ (dX) } V, = [d;-k^f V, , 

(2.33) 

= (l/(/yc)) {-Re{/:* (F,,-c)v,+ (i/^v,) v.} + [£/^^-A:']vJ. 

(2.34) 

Boundary conditions for both Couette flow as well as Poiseuille flow 

quantities are given below: 

At fluid solid interface (at z = 1): 

in term of perturbed 

v,=0. 

(2.35) 

v,=0. 

(2.36) 


At fluid-deformable wall interface (at z= 0): 
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For Coiiette flow -zand for Poiseuille flowv^ =z (1 — z) , the second term of left hand 

side of the tangential velocity condition given below is obtained by linearizing the x- 
component interface velocity about the unperturbed interface at z = 0: 


V, = -ikcu^. 

(2.37) 

v\+ li, = -ikcu^, 

(2.38) 

4 = <?.xz > 

(2.39) 

II 

(2.40) 

ikv^. 



-P/+ 2 

^ 2 {^-ri^ikc)d^u^, 


<?zz = (-p-^r ikc)[dji^ + iku^]. 

At rigid solid-deformable wall interface {z =-H): 

u,= ^, (2.43) 

(2.44) 

We solve above equations using a combination of asymptotic and numerical methods. 
First we would like to do the asymptotic analysis for high Reynolds number {Re) and 
then we will continue asymptotic results to moderate values of Re using numerical 
technique. 
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2.3.2 Numerical method for solving characteristic equations 

There are many methods by which Orr-Sommerfeld equation can be solved but the 
stiffness of the equation at high Reynolds number requires extra care in the solution. The 
main aim of the numerical solutions for a given basic flow are (1) to find the curve for 
marginal stability (c,. = 0) and the curve of constant growth rate, (2) to find the eigenvalue 

spectrum for a given pair of positive values of c. and Re, and (3) to calculate the 
associated eigenfunctions. 

There are three accurate methods to solve Orr-Sommerfeld equation those are 
Spectral expansion, i.e. expansion of solution as a linear combination of a complete set of 
orthogonal functions, use of finite differences to integrate the Orr-Sommerfeld equation 
as initial value problem whereby integration is from one boundary towards the other. We 
use the principle of expansion of the solution as a linear combination of orthogonal 
solutions to the governing stability equation which is simple to use as well as accurate. 

In our case, we have fourth order ordinary differential equation for both fluid and 
deformable wall as given in equations 2.33 and 2.35. Hence to know the complete 
solution of fourth order ODE we must know functional value and their first, second and 
3'^^' derivative. This means we need to choose four orthogonal functions and hence four 
constants to write general combination of those four functions. 


Qf — A-\ 01+ A 2 02 
Fluid 


Interface 


/ / /,' OfefQf'mable Wall,-"' ° 

/0„/=Bi/0it'B2 02 / / / / 


Solid 


Fig. 2.4 Schematic diagram for numerical calculation 
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However the boundary conditions at both the rigid walls allow us to take only two 
orthogonal functions satisfying two given boundary conditions 2.37 and 2.38 at upper 
rigid wall and conditions 2.43 and 2.44 at lower rigid wall. Thus we need only two 
constants to write solution set for velocity of fluid and two constants for the deformation 
field. After that we used fourth order adaptive Runge Kutta as an integrator and marching 
from the rigid surface to interface of the fluid and deformable wall. We have four 
interface conditions 2.39, 2.40, 2.41 and 2.42, and after applying interface conditions we 
are able to get four linear homogenous equations with four constants. As we know for 
getting non-trivial solutions (i. e. non-zero value of A|, A 2 , Bi and B 2 ) of homogenous 
equations, determinant of the coefficient matrix must be zero. Hence by setting 
determinant equal to zero, we obtained the characteristic equation. This was solved 
numerically using Newton-Raphson, iteration technique. We used the programming 
language ‘C’ for developing code. To compute numerical we must have to know 
sufficiently close initial guesses and for that we used asymptotic results for high Re (V. 
Shankar personal communication). 

Mathematically the above procedure can be represented as follows: 

According to the mentioned procedure two independent orthogonal solutions for the 
velocity field of fluid consistent with the BC at z = 1, can be written as: 


0 

II 

0 

II 

di^^=\ 

II 

0 

v:=o 

0 

11 

II 

0 



The velocity field is a linear combination of these two solutions. We can write solutions 
set in combined form as Q = d^v^, d^v^, dlv,) where 0i={O, 0, 1, 0} and 

02= {0,0,0, 1}. Similarly for deformable wall the displacement fields at z = -H, 
consistent with the BC at z= -H: 



= 0 

d^u^ = 0 

11 

c/jM. = 0 


= 0 

i! 

0 

= 0 
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Here {(), 0, 1, 0}and 02“{O, 0,0, 1} and by using matching condition at the 
interface of fluid and deforinable wall, we will be able to get the characteristic equation, 
which has all the information about the stability of the system of interest. 

23t3 Validation of numerical method 

lo check the validity of numerical solution we compare our numerical results with the 
result of [1999] Kumaran for low Reynolds number flow past a deformable wall and 
Shankar and Kumaran ’s wall mode numerical results and asymptotic results. Although 
our code matches with both the results, I have shown comparison only with Shankar and 
Kumaran's wall mode results. 



Re 


Fig. 2.5 r vs. Re in Couette flow past a deformable wall: Data from numerical solution 
for //= 2, k = 0.5, 0 for first two modes. The two dotted lines are asymptotic results. 
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Re 

Fig. 2.6 r vs. Re in Couette flow past a deformable wall: Data from numerical solution 
for H- A: = r]r- 0 for first two modes. The two dotted lines are asymptotic results. 

2.4 Results and Discussions 

Z4.1 Couette flow past a deformable wall 

The transition strain rate at which the perturbations become neutrally stable (c,- = 0) 
depends on the wave number k, the wall thickness H, the viscosities ratio rjrand Reynolds 
number Re. In this section, the behavior of the transition stain rate at 17 ^ = 0 is analyzed. 
We already mentioned in introduction that there are several neutral modes in wall mode 
asymptotic analysis. Hence the corresponding numerical solution also has those modes. 
Out of that, first mode is viscous mode same as presented by Srivatsan and Kumaran, 
(1997) with slope -1/3 at high Reynolds number; it also gives the results presented by 
Kumaran (1994) at low Reynolds number. But when we continue the second wall mode, 
it first has slope -1/3, which is expected since it is continuation of wall mode. But at 
lower Reynolds number its turns back with slope -1 and it continues to intermediate as 
well as high Reynolds number. This can be observed clearly in Fig. 2.7 and Fig. 2.8. 
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Gamma 



Fig.7 r vs. Re for Couette flow past a deformable wall with //= 1, 

1e+08 I - 1 - — - r-- 


k= 1, rir = 0. 

model 
mode2 



1e-06 i ^ ' — ■■■' ^ " ^ ^ ' 

0.0001 0.001 0.01 0.1 ^ 10 100 1000 10000 100000 1e+06 


Re 


Fig-2.8 r v$. Re in case of Couette flow past a deformable wall for //= 2, k = 0.5, A]r= 0 
Data from full numerical solution. 


Gamma 


The transition strain rate Ft is shown for //= 1, A: = 1, = 0 and for 2, yt = 0.5, 

nr=0 as a function of Reynolds number in Fig. 2.7 and 2.8 respectively. Fig. 2.7 and 2.8 
reveal that most rmstable mode corresponds to scaling F ~ Re'^ , which does not have 
same scaling as wall mode. Thus, it must be a different mode since wall mode has scaling 

- 1/3 • • 

F ~ Re . Later in this section we will inspect these modes extensively. It can also be 
observed from Fig. 2.7 and 2.8 that mode-2 has two different curves. 



Re 

Fig.2.9 F vs. Re in case of Couette flow past a deformable wall for 77= 1, k = 1, A7r= 0 
from full numerical solution corresponding to second asymptotic wall mode. 

As we have shown in Fig. 2.9 wall mode-2 provides two different curves as 
compared to all other modes, which have only one. Thus we would like to explain the 
reason for obtaining two curves for this particular mode. We already had mentioned that 
when we continued wall mode-2 to lower Reynolds number then it turns back and 
continues to high Reynolds number unlike wall mode-3 and wall mode-4 in Fig. 2.7 and 
Fig. 2.10. 
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Fig.2.10 r VS. Re in case of Couette flow past a deformable wall for k = 1, /7r = 0 
from full numerical solution corresponding to second asymptotic wall mode. 

We further analyzed the specific region to check the domain of instability. We 
found that the region above upper curve is unstable^ regime within the U-shape curve is 
stable and region below the lower cure of Fig. 2.10 is unstable. But we know that wall 
modes are always stable in the case of flow past a rigid wall. However, the curve in Fig. 
2.10 contradicts this expectation since at vanishing F deformable wall behaves like rigid 
wall. This reveals the possibility of an extra curve and therefore we tried to explore the 
whole lower region of Fig. 2.10 by lowering F. We found that after certain values of F 
imaginary part of growth velocity is changing from positive to negative. This indicates 
the possibility of another neutral curve which appears as the lowest curves which was 
there in Fig. 2.9 and hence we got stable region as given in Fig. 2.11. We already had 
mentioned about slope of every part of the figure and also discussed that new part must 
be a different mode than wall mode. 
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Gamma 


1e+08 



Re 

Fig.2.1 1 r vs. Re in case of plane Couette flow past a deformable wall forH= 1, 
k= 1 , / 7 r = 0 for second wall mode and their extension. 

When one tries to find the instability experimentally then one will able to get only 

the most unstable mode, which is not the wall mode. So it is instructive to examine the 

eigenfunction of the neutral modes obtained from numerical solution in order to find the 

validity for wall mode and to find causes for new modes. We find the wall mode has 

same scaling as my Guide Shankar had used for asymptotic analysis and also Shankar 

and Kumaran [2002] presented. But when we further analyzed the eigenfunction for 

lower modes then we found that it has different scaling than wall mode. As mentioned 

above F scales as Re’' as compared to wall mode where F scales to Re'^^ and the real part 

1/3 

of growth velocity remains constant as compared to wall mode where it varies as Re' 
and further we can get more idea about the nature of new mode by knowing scaling for 
all the dynamical quantities with Re as well as velocity distribution as a function of 
position. Thus we plot velocities and their derivatives as function of Reynolds number on 
log-log scale as presented in Fig. 2.12-2. 17. 
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Amplitude of velocity functions and derivative 


Proni above mentioned figures we are able to obtain scaling as shown in table 1 
for the quantities related to fluids and in table 2 for the deformable wall. 



Vz 

Vx 

dzVx 

r 

Pr 

Cr 

Wall mode 

Re" 


Re"'" 

■ Re-'^" 

Re' 


Lower mode 

Re" 

R? 

Re''" 

Re' 

Re' 

-R? 

The lowest 

mode 

Re" 

Re" 

Re''" 

Re' 

Re' 

“r? 


Table! . Scaling of different quantities in fluid as a function of Reynolds number 
for H = 1 , k = 1 , //r = O.Obtained from our numerical code. 



Re 


Fig. 2.12 Variation of the absolute value of velocity field and their derivative in 
the fluid with Reynolds number for wall mode with //= I, A: = I, /7r= 0. The two dotted 
lines are reference line with slopes 1/3 and 2/3. 
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Amplitudes of velodfy functions and derivative Amplitude of wall displacement fields and derivative 



I' ig. 2. 1 3 V ariation of the absolute value of displacement field and their derivative 
in the fluid with Reynolds number for wall mode with H = 1, k = 1, 17 ^ = 0. The three 


dotted lines are reference line with slopes 1/3, 1/2 and 1/3. 
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Fig. 2.14 Variation of the absolute value of velocity field and their derivative in 


the fluid with Reynolds number for lower curve with H- 1, /c - 1, f 7 r - 0. 
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Amplitudes of velocity functions and derivative 



Re 

f'ig. 2.15 Variation of the absolute value of velocity field and their derivative in 
the fluid with Reynolds number for the lowest curve with H= \,k= \, qr= Dotted line 



F‘ig. 2. 1 6 Variation of the absolute value of displacement field and their derivative 


in tlie fluid with Reynolds number for lower curve with //= I , A- = 1 , /7r = 0 
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Amplitude of displacement field and 



Re 

Mg. 2. 1 7 Variation of the absolute value of displacement field and their derivative 
in the fluid with Reynolds number for the lowest curve with H= k= 1, 17^= 0. The 
three dotted lines in Fig 2. 1 6 and 2. 1 7are reference line with slope 0. 




Ux 

dzUx 

Ft 

Pg 

Cr 

Wall mode 

Re^ 

Re^^^ 


Re-'^' 

~R? 

Re'^^ 

Lower mode 

Re" 

R? 

Re" 

Re-' 

Re' 

Re" 

The lowest 

mode 

Re" 

Re" 

Re^ 

Re' 

Re' 

Re" 


Table2. Scaling of different parameters in deformable wall vs Reynolds number for H = 
l,k= 1, 0. 

We can see from tables the scaling of parameters tells us that lower curve 
corresponding to mode-2 has exactly same scaling as the lowest curve. This reveals that 
these two inodes have same properties. This scaling is similar to inviscid mode at high 
Reynolds number mentioned by Shankar and Kumaran [2002]. Further results can be 
confirmed by the analysis of distribution of velocity field and displacement field 
(eigenfunctions) as a function of position. Other than this boundary layer thickness vs. 
Reynolds number graph also reveals that the lower and the lowest modes are inviscid 
modes. 
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Amplitude of boundary layer 


1 



z 

Fig. 2.18 Variation of thickness of wall layer (Boundary layer) with Reynolds 
number for plane Couelte (low with H= \,k= \,nr=Q. The three dotted lines have slope 
-1/3, -1/2 and -1/2. 

Figure 2.18 shows the variation of the thickness of the wall layer S with the 
Reynolds number. The wall layer thickness may be estimated by computing the ratio 
v/c4v, at z == 0 from the full numerical solution. As Fig. 2.18 shows, the wall layer 
thickness S decreases as Re'^'^ for upper curve and as Re^ for other two curves in the 
limit of large Re. Upper curve results matches with the asymptotic analysis for the wall 
mode case but other two have different scaling which again tells us that these modes are 
inviscid modes. This shows that the numerically observed upper curves are indeed wall 
mode and lower curves are inviscid modes. For further conformation we can calculate 
eigenfunctions. 

7'o calculate eigenfunctions we must know the con.stants multiplying the linearly 
independent numerical solutions, which are mentioned in Fig 2.4. But the equations we 
obtain after substituting in the boundary conditions are homogenous and linearly 
dependent. So it is not possible to find a unique solution, as the determinant of coefficient 


28 



|vz| 


matrix is zero for linearly dependent equations. Thus we change one of the equations 

with (Vz)j=o — (1 + 0/^/2 and solve for constants. By this, we are fixing the magnitude of 

fluid velocity. Then we are able to get the distribution of eigenfunctions. We got same 
behavior for wall mode as presented by Shankar and Kumaran (2002). 



0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 


z 

Fig.2.19 The amplitude of velocity |Vz| vs. z in case of Couette flow fox H=2,k = 
0.5, nr =0 at Reynolds number 3000 corresponding to all the curve of mode 2. 

Figures 1.19 and 2.20 show the variation |vzl and |vx| as a function of z in the fluid, 
and egienfunction are calculated subjected to normalized condition mentioned above 
paragraph. These figures show that corresponding to upper curve behavior is same what 
Shankar and Kumaran (2002) had presented in their work. But for lower cures velocities 
are uniformaly-distributed through out the channel and there is large variation near the 
wall. This also confirms that these new modes are inviscid mode. 
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Fig.2.2() The ainplitucle of velocity |Vx| vs. z in case of Couette flow for Tf = 2, A: = 0.5, nr 
=0 at Reynolds number 3000 corresponding to all the curve of mode 2. 

Further we can analyze the flow for high Reynolds number limit using the 
obtained scaling for the new mode which is claimed to be an inviscid mode. When we 
perform asymptotic analysis by using new scaling, we got equations as follows which 
further conform that obtained new modes are inviscid mode. 


ik + {ci^v^)v, - ~ik Re"' , 

ik (v, -c)y. == - Re ' . 

- k\-^u, = -ik Re"' [dyk^]u^^ 

- kVu, = - Re ■' r/,/T, +— [d^ -k^]u , . 

*• ^ J"* • * 


(2.45) 

(2.46) 

(2.47) 

(2.48) 
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Gamm 


2.4.2 Results and discussion for the Poiseuille flow past a 
deformable wall 

We already had discussed about the results of plane Couette flow and their 
numerical validity. Theie we confirm that both viscous as well as inviscid modes are 
present in How past a deformable wall. For plane Poiseuille flow we just have to change 
base velocity profile from plane Couette to plane Poiseuille in Orr-Sommerfeld equation 
in a way mentioned in section 2.3.2. Besides that there must be an extra term in the case 
of plane Poiseuille flow in Orr-Sommerfeld equation since in this case second derivative 
of base velocity profile is non zero as compared to plane Couette flow. Other than above 
changes, all the terms remain same for both the cases since we required only first 
derivative at ^ 0 which is same for both the cases. The asymptotic results also agree for 

high Reynolds number, which is same for both the cases. 



Re 


Fig.2.21 Gamma vs. Re in case of Poiseuille flow past a deformable wall for // = 
\,k~ \ ,rir~(): Data from full numerical solution. 
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Amplitudes of displacement 


In figuie 2.21 we had plotted F as a function of Re and the slope of the line at 
high Reynolds number reveal the existence of wall modes as well as inviscid mode as 
found in Couette flow. But for intermediate Reynolds number it has very different 
behavioi than Couette flow especially for wall mode 1 and 2. As compared to plane 
Couette flow in Poiseuille flow first wall mode does not capture low Reynolds number 
mode and when we extended second wall mode to intermediate as well as low Reynolds 
number then instead of converging to inviscid mode it continues to low Reynolds number 
and shows same behavior as other higher modes. When we again tried to understand the 
behavior of low Reynolds number mode at intermediate and high Reynolds number then 
it continues to the invi.scid mode. To get final confirmation about these modes we 
determine scaling of various quantities as a function of Reynolds number and presented 
in tables .1 and 4. from these tables we confirm that these are the same modes what we 
had anticipated above. 



Re 

Fig. 2.22 Variation of the absolute value of velocity field and their derivative in 
the fluid with Reynolds number for the wall mode Poiseuille flow with FI= 1, A:= 1, rjr — 
0. The two dotted lines in are reference line with slope 1/3 and 2/3. 
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Amplitude of veiocity 


100000 



Re 

I'ig. 2.23 Variation of the absolute value of displacement field and their derivative 
in the Huid vs. Reynolds number for wall mode Poiseuille flow with H= \,k=\, r]r- 0. 



Fig. 2.24 Variation of the absolute value of velocity field and their derivative in 
the fluid with Reynolds number for inviscid mode Poiseuille flow with H= rjf = 

0. 
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Amplitudes of displacement 



Re 

Fig. 2.25 Variation of the absolute value of displacement field and their 
derivatives vs. Reynolds number for the inviscid mode in Poiseuille flow with H= \,k = 

l,'7r=0. 



0.0001 

1000 10000 Re 100000 1e+06 


Fig. 2.26 Variation of thickness of wall layer (Boundary layer) with Reynolds 
number for Poiseuille flow with //= 1, /r = 1 , i 7 r= 0. The two dotted lines have slope -1/3 
and -1/2. 
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Tables. Scaling of different parameters in the fluid past a deformable wall with 
Poiseuille base velocity profile flow vs Reynolds number for H=\,k= I, r]r=0. 


Wall mode 



Table4. Scaling of different parameters in the deformable wall with Poiseuille base 
velocity profile flow vs Reynolds number fox H= \,k = rir = 


1 



0.4 0.6 

Z 


Fig.2.27 The amplitude of velocity lVz| vs. z in case of Poiseuille flow for 1, * = 1, 
/ 7 r= 0 at Reynolds number 1000. 











2.5 Conclusion and future work 

1 he stability of wall modes and inviscid modes in the fluid flow past a deformable wall is 
analyzed using a combination of asymptotic and numerical method. In order to 
demonstrate the instability, we used the case of Couette flow as well as Poiseuille flow 
past a deformable wall. We have shown that wall mode instability predicted by the 
asymptotic analysis is also captured by numerical methods for both Couette flow as well 
as Poiseuille flow past a deformable wall. Further we have shown that there exists a mode 
in Couette flow which has different scaling than that of wall modes and this mode is also 
present in the Poiseuille flow. But there is a difference between in these two cases. In 
Couette flow one of the wall modes captures the inviscid mode when we extend wall 
mode solution numerically. But in Poiseuille flow instead of wall mode, zero Reynolds 
number modes capture the inviscid mode. 

To examine the type of modes, we calculated the scaling of velocities and their 
derivatives with Reynolds number for the case of Couette flows and we find that new 
modes are indeed inviscid modes. Further we find the boundary layer thickness as a 
function of Reynolds number and this also indicates that the new mode is an inviscid 
mode. Other than these, when vye find velocity distribution in the channel and it shows 
that there is uniform distribution in the channel as compared to wall modes where the 
velocity field is confined near the wall which confirms that this new mode must be an 
inviscid mode. But for the case of Poiseuille flow, we got the same scaling as Couette 
flow but velocity field distributions have a different nature instead of expected critical 
layer for invi.scid modes. This prediction of the present study needs further investigation. 
Other than this, our numerical results predicted the inviscid modes for both the plane 
Couette as well as plane Poiseuille flows. The prediction of inviscid modes for plane 
Couette flow in the present study is in stark contrast with the asymptotic prediction of 
Kumaran [1994b], who concluded that the inviscid modes are stable in the case of plane 
Couette flow. 
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CHAPTER 3 


Linear stability analysis of two-layer plane Couette flow 
past a deformable solid layer 

3.0 Introduction 

The study of stability of two-layer fluid flow past a deformable wall is important because 
of its use in multi-layer extrusion, where an accurate understanding of stable and unstable 
processing conditions can help in preventing interfacial instabilities, which is important 
to achieve desired product quality. In this chapter, we intend to study the effect of a 
deformable solid layer on the stability of two-layer plane Couette flow in a channel. We 
are motivated to study this problem by considering the possibility of coupling between 
two interfacial waves of the fluid-fluid & fluid-deformable solid wall interfaces. First we 
analyzed asymptotically for long wave and zero Reynolds number limits and then 
continued the asymptotic results numerically to intermediate values. Before our study on 
stability of two-layer fluid flow, first significant study had been done by Yih [1967]. Yih 
carried out an asymptotic analysis of the interfacial mode for long waves and showed that 
a jump in viscosity across the interface can lead to long wave instability. Moreover Yih 
find that long wave instability is present when more viscous fluid is thinner otherwise 
long wave is stable. 

Subsequently Hooper and Boyd (1983) solved the linear stability 
problem for two-dimensional plane Couette flow using a short wave asymptotic analysis. 
They showed that this flow is unstable for short waves as well. The asymptotic 
expression for the growth rate of the interfacial mode for short wave is essentially 
consists of three terms. Out of that dominant term associated with surface tension, which 
is always stabilizing. The next term depends on gravity and the ratio of the densities and 
third term depends on ratio of viscosities: this term is always destabilizing if the density 
are equal. 

In present study, we first analyze for the problem for fluids with equal density 
then we will examine the effect of density difference on the stability of fluid-fluid 
interface. Hence out of three-mentioned terms, for short waves as shown by Hooper and 
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Boyd, only first and third are valid for matched fluid densities. The case of equal 
densities of the fluid, which is om case there exist some rules of thumb given in text book 
Fundamental of two-layer Dynamics by Joseph & Renardy, which tells that we can 
stabilize long wave by placing the less viscous fluid in a thin layer and stabilize the short 
waves by assuming sufficiently high surface tension. In this chapter, we illustrate the 
following effects; 

1) Two-layer flow past a rigid wall is unstable and it is possible to suppress the 
instability by making the solid layer sufficiently deformable. 

2) Two layer flow past a rigid wall is stable and it is possible to destabilize the flow 
with deformable wall and, 

3) We will examine the effect of various parameters on the stability of both the 

interfacial modes. 

3.1 Problem formulation 

Our system consists of two fluid layers of thickness P R and (1-P) if respectively 
and the base is supported by a deformable solid of thickness H R and upper layer is 
bounded by a rigid wall. The fluids are assumed to be incompressible, Newtonian fluids 
top rigid wall having shearing tangential velocity equal to U*. Governing equations for 
both fluids and deformable solid for mass and momentum transfer in the bulk phases are 
the continuity equations and the Navier-Stokes equations. Schematic diagram for above 
situation is shown in Fig. 1. 



Fig 3.1 Schematic diagram of two-layer plane Couette flow past a deformable wall 
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3,2 Governing equations 

3.21 Governing equations in non-Dimensional form 

The basic hydrodynamic equations and non-dimensionaiisation parameter remains same 
as in chapter 2 except, for dimensionalisation of viscosities and densities, those are 
dimensionalised by viscosity of fluid- 1 and deformable wall density respectively as 
follows: 

= p‘‘=p‘‘lp\ 

y = 

Here a = 1, 2. p*, y* and y are dimensional density of deformable solid, dimensional and 
non-dimensional surface tension between fluid-fluid interface respectively, p*" are the 
dimensional density of two fluids, p^, //, and are the viscosities of deformable solid, 

fluid-1 and fluid-2 respectively. 

Thus here we can start directly with non-dimensional form of equations. 

Mass balance for fluids: 

5,vr=0, (3.1) 

Mass balance for deformable wall: 


d^u, = 0 , 

Momentum balance for fluid and deformable wall can be written as: 


(3.2) 


Re p'^id, + v;5,)vr = -5,;?; ■ 

Re(aV^, ) = -Pg^,j+ 


(3.3) 

(3.4) 


Here ar=l, 2 represents the fluid index. Here repeated indices are summed over 
and single index means component of a vector. Here v, u, Pf- and pg represents velocity 
of fluid, deformation in viscoelastic wall, pressure inside fluid and pressure inside 
deformable wall respectively. 
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The no-slip condition at rigid solid and fluid interface is (at z =1); 

vf = 0, 

The impermeability condition at the solid surface is (at z = 1): 

vf) = 0, 

Continuity of velocity component and stress component can be written as: 






( 3 . 5 ) 

( 3 . 6 ) 

( 3 . 7 ) 

( 3 . 8 ) 

( 3 . 9 ) 


n,Tpn,= n,Tpn,, ( 3 . 10 ) 

The tangential and normal velocity balance at the fluid-defoiTnable wall interface: 


= d.ii. 


vf = , 


The shear stress balance at fluid- deformable wall interface; 

/! = t cy n 

The normal stress balance at fluid-deformable wall interface: 

The no slip-condition at deformable wall-rigid solid surface is (atz =-// ); 

“x =9. 

The impermeability condition at deformable wall-solid surface is (atz =^/7 ): 


=0 , 


( 3 . 11 ) 

( 3 . 12 ) 

( 3 . 13 ) 

( 3 . 14 ) 

( 3 . 15 ) 

( 3 . 16 ) 


Non-dimensional stresses in liquid and deformable wall can be given respectively as 
follows: 
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t;; =-p;<5, +/.ro.v; + dX), 

= ~^iPK + (p + l.d,){diU, + d,u.). 

Kinematical equation for fluid-fluid interface can be written as; 
d,h + dji= . 

Where h is the height of the fluid-fluid interface from the x-axis. 

3,Z2Base state calculation for both Poiseuille and Couette flow: 

After solving two-dimensional Navier-Stok.es equation in non-dimensional form for 
above situation we obtained the base velocity profile for fluids for the case of equal 
density are as follows; 

V,^ =(z-^-^(Ur-V)/0+P(lJr-V); 

V,^ =Z^Jr/(^^^3(lJr- V); 

V,, =0; 

Vj. = 0; 

Base state deformation for gel in non-dimensional form is; 

3.3 Linear stability analyses 

3.3t1 General equation for linear stability analysis 

Base state for Couette flow is given in section 3.2.2. We now add perturbation in base 
profile to find out hydrodynamic equations and boundaries condition in terms of 
perturbed quantities. 

Thus perturbed quantities can be represented as follows; 

-Tf -Po^P,- = [vj,r^,u,p",p^,s] exp{ik {x-ct)} 
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Where, c and k are the growth velocity and the wave number of the instability 
respectively and base state condition is given in section 3.2.2. '’s' represent perturbation 
amplitude in the height of fluid-fluid interface. 



Fig 3.2 Schematic diagram of two layers plane Couette flow past a deformable wall in 
perturbed state 

Applying above perturbation to equation 3.1 to 3.16 and equating linear term on both the 


sides, we obtain: 

/A:v;=0, (3.17) 

Ktp^ik (v; -c)F; + = -ikp- +I.C [d] -k^K . (3-18) 

Kep^ik (F;-c)Fn = [dl-k^K ■ (3-19) 

For deformable wall continuity equation, x and z component momentum balance for 
perturbed states are respectively: 

d,u,+ iku^ = 0, (3.20) 

-Rtec\=-ikp^+i^-rjpkc)[dl-k^]u^, (3.21) 

-Re + (-p-??,. /^c)[(7.^ . (3.22) 
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After combining above equations 3.17 to 3.19 for fluid and 3.20 to 3.22 for deformable 
wall we get a standard equation, which is known as Orr- Sommerfeld equation for fluid 
and deformable wall respectively. 

For deformable wall: 


-Rq k^c\dl-k']u^= (3.23) 

pg = (1 /(/ A:)) { Re + (^ - 7,. ^ ^ } . (3 .24) 

For fluid: 

Re p-‘it{(v;-c) [dl-e]+ (d,V;) }vf =< [dl-efv‘, (3.25) 

p-, = (ll(ik)){-Ktp‘{ik (vf-c) 5“ + (rf.v;) v“}. (3.26) 

Boundary conditions in term of perturbed quantities can be written as: 

At fluid solid interface (at z = 1): 

= 0 , (3.27) 

vf^ = 0, (3.28) 

At fluid-fluid interface (at z=p): 


The second term of left hand side of the tangential velocity condition given below is 
obtained by linearizing x-component interface velocity about the unperturbed interface at 


z z 

(3.29) 


(3.30) 

II 

(3.31) 

fS'’ = fS>-rk‘s 

(3.32) 


At fluid-deformable wall interface (at z= 0): 
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The second term of left hand side of the tangential velocity condition given below is 
obtained by linearizing ^-component interface velocity about the unperturbed interface at 
z = 0: 


-ikcti ^ , 

+ (^z'^r'^)z=o “z = -ikcu^, 

^i”=‘?zz, 

At rigid solid-deformable wall interface (z =-//) 

M, = 0, 

u^= 0 . 

Here, 

f:=-Pf+2/jA^% 

A = -Pa 2 . (p-77, ikc)d^u^ , 

ikc)[dM, +iku^], 

-ikes + vl'\s = Vj'' . 

We can solve above equation using a combination of asymptotic and numerical methods. 
First we would like to solve for some asymptotic cases for long wave and zero Reynolds 
number limit then we will continue asymptotic results to moderate value using numerical 
technique. 


(3.33) 

(3.34) 

(3.35) 

(3.36) 

(3.37) 

(3.38) 
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3.3.2 Numerical method for solving characteristic equations; 

In our case we have fourth order ordinary differential equation for both fluids and 
the deformable wall as given in equations 3.23 and 3.25. Hence to know the complete 
solution of fourth order ODE we must have to know the function’s value and their first, 
second and 3*^*^ derivatives. This means we need to choose four orthogonal functions and 
hence four constants to write general combination of those four functions. 


Solid wall 


©2 — A 2 02 

Fluid 2 


©1 =/\3 03+/\404 +/\5 03 + Aa 06| 

Fluid 1 


Fluid-deformable wall Interface 


/ / / / / / Deformable wad / / 

t t / / / / // / r i 

/'0i=,Si'(^i+.e2^2/ / / / / / / / 


/ 0 


Solid wall 


Fig. 3.3 Schematic diagram for numerical calculation 

The boundary conditions at both the wall allow us to take only two orthogonal functions 
which satisfying two given boundary conditions 3.27 and 3.28 at upper rigid wall and 
conditions 3.37 and 3.38 at lower rigid wall. Thus we need only two constants to write 
the solution set for velocity of fluid-2 as well as deformation field but we must need four 
functions to represent complete solution for fluid- 1. After that we used fourth order 
adaptive Runge Kutta as an integrator and marching from the rigid surface to interface of 
the fluid and gel. We have eight interface conditions 3.29 to 3.36, after applying these 
interface conditions we are able to get eight linear homogenous equations with eight 
constants. As we know for getting non trivial solutions (i. e. non-zero value of A|, A2, A3, 
A4, As, Ae, Bi and B2) of homogenous equations, determinant of coefficient matrix must 
be zero. Hence by setting determinant equal to zero, we obtained the characteristic 
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equation. This was solved numerically using Newton-Raphson, iteration technique. We 
used programming language ‘C’ for developing code. 

Mathematically above statement can be represented as: 

According to mentioned procedure two independent orthogonal solutions for the velocity 
field of fluid-2 can be written as: 


v,=0 

d^v,^0 

dlv^ = \ 

r/^'v,=0 

v,=0 

d^v^=0 

o 

11 

dlv^ — \ . 

u 

11 

d^v^=0 

dlv,={) 

d]v.=Q 

o 

II 

N 

dy^ = 1 

dlv^ = Q 

dlv=0 

O 

II 

N 

d^v=0 

dlv, = \ 

y.\ 

II 

o 

v.=0 

dv^ =0 

dlv^=Q 

d]%=\ 


The velocity field is a linear combination of these four solutions. We can write solutions 
set in combined form as 0= {v^, where 0i = {0, 0, 1, 0} , 

02= {0, 0, 0, 1} , 03= {1, 0, 0, 0} , 04= {0, 1, 0, 0} , 05= {0, 0, 1, 0} and 06= {0, 0, 0, 1} . 
Similarly for deformable wall the displacement fields are: 

= 0 = 0 d^ii^ = 1 d^u^ = 0 

u, = 0 d/i^ = 0 d^u^ = 0 r/Jw. = 1 . 

Here 0i = {0, 0, 1, 0}and 02={O, 0, 0, 1} and by using matching condition at the 
interface of fluid and deformable wall. We will able to get characteristic equation which 
has all the information about the stability of the system. 
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3.4 Results and discussions 

3 , 4.1 Numerical validation and results for Rigid channel 

To check the validity of numerical solution we need to put the value of = 1 and surface 

tension zero in current numerical code. By putting these values two-fluid problem 
becomes single layer flow past a deformable wall, which is rightly able to capture all the 
results presented in Chapter-2. First it is instructive to briefly recapitulate the results for 
the stability of two-layer Couette flow past in a rigid channel and this has been the 
subject of earlier study of Yih [1967] and Hooper and Boyd [1983]. Yih carried out 
asymptotic analysis for long wave using a formal power series in the wave number k « 1 
for k Re « 1, where Re is Reynolds number and k is wave number. Yih show that jump 
in viscosity across the interface can lead to long wave instability which can be observed 
from Fig. 3.4 and 3.6 which is computed numerically taking asymptotic results as a initial 
guess. 



k 


Fig. 3.4 Fluid-fluid interface instability in rigid channel (r= 0): variation of the 
imaginary parts of the wave speed c with the wave number k for different values of 
thickness ratio p, interfacial tension = 0, viscosity ratio Mr ~ 0-5 and Reynolds number Re 




0.001 0.01 0.1 1 10 

k 


Fig. 3.5 Fluid-fluid interface instability in rigid channel (F = 0): variation of the real parts 
of the wave speed c with the wave number k for different values of thickness ratio p, 
interfacial tension = 0, viscosity ratio 0.5 and Reynolds number Re = 1. 

Further Hooper and Boyd [1983] performed an asymptotic analysis for short 
waves and showed that short waves are unstable and surface tension has ability to 
stabilize the instability associated with short waves, which can be easily observed by 

seeing the difference between Fig 3.4 and Fig 3.6 is plot between C. vs. k with surface 

tension zero thus it is unstable for all ranges of wave number similar to what Yih and 
Hooper & Boyd found in their work. But in Fig 3.6 flow is unstable for small wave 
number and gets stabilized for long wave number due to presence of surface tension as 
mentioned by Hooper and Boyd [1983]. Although these authors talk of short and long 
wave number, the most unstable mode has finite wave numbers as observed in Fig. 3.4 
and 3.6. Thus we must need full a numerical solution to analyze the stability for whole 
range of wave number including intermediate wave numbers. 
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0.0005 


P =0.4&H=0.5 
P =0.2&H = 0.5 



Fig. 3.6 Fluid-fluid interface instability in rigid channel (r= 0): variation of with k for 
different value of P, interfacial tension = 0.01, //,, = 0.5 and /?c = 1. 



k 

Fig. 3.7 Fluid-fluid interface instability in rigid channel (r = 0): variation of c,with k for 
different value of P, interfacial tension = 0.01, 0.5 and Re= \ . 

VTRcfttT jPTJrr* 

wf%*roA...IiSl8i„L 
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3.4.2 Asymptotic results 

Before proceeding further we want to summarize the asymptotic results (V. Shankar, 
personal communication) showing the effect of F on the stability of our system of flow 
for the case of long waves where analytical results are available. Analytical results for the 
long waves are obtained as follows: 

c=c, + ikc^, (a) 

Where is a real number does not depends on F &. Re and c, is also a real 
quantity know as first correction factor which is responsible for the stability and 
instability of flow is a function of F as well as Re ( whenc^X) flow is unstable and c, <0 
flow is stable). We first had presented an analytical results for interfacial tension = 0.1, 
= 0.5 and H - 0.5 as in equation (b) for the case when two-layer flow is unstable in 
rigid channel and deformable wall has stabilizing effect. This can be observed by 
equation (b) that F = /< U* /G R always be non-negative hence it has always stabilizing 
effect other than the case of rigid wall for which F is zero (elasticity of the solid layer is 
very high). 

c= 0.767023 + / (-0.00527865 *F + 0.000 124662* Re) (b) 

Again we would like to discuss an analytical result for the case when two-layer 
flow in a rigid channel is stable and deformable wall has destabilizing effect. This is 
possible for the case of thin layer effect (i.e. thinner fluid has higher viscosity). For 
obtaining analytical results for this case we used the variable quantities P = 0.8, 
interfacial tension = 0.0, //,.= 0.25 and //= 2.0 and will able to get expression given in 

equation (c). From equation (c) it is obvious that for this case F has a destabilizing effect 
and flow remains stable for the particular value of F know as transitional value but if we 
will further increase the value of F than our system becomes unstable as suggested by 
equation (c). In proceeding parts of the thesis we intended to examine the effect of 
deformable wall on the stability of flow for moderate and high values of wave numbers. 
c=0.720214+/ (0.44450 * F- 0.000155329* Re) (c) 


52 



3.4.3 Fluid flow past a rigid wall is unstable and deformable wall has 
stabilizing effect (/^^ > 1 & p > 1/2 or < 1 & p < 1/2) 

In this section we intended to analyze the arrangement of the two-layer flows, which is 
unstable in rigid channel and deformable wall has stabilizing effect. From Fig 3.8 it is 
clear that it is possible to stabilize low wave number instability by changing rigid wall to 
deformable wall as presented in equation (b). This can be explained by considering the 
interference between two interfacial perturbations namely fluid-fluid and fluid 
deformable wall. Further it is clear that it is not possible to stabilize short wave length 
instability irrespective of F value. This is because; there is no possibility of coupling 
between two interfacial perturbations for short wave lengths (high wave number 
perturbations) as they are confined near the interface only. Other thing we can observe 
from the same Fig. is first stability increases as elasticity of deformable wall (T value) 
increases but after certain limiting value of T , if we further increase the value of T flow 
becomes unstable. But we know from the Hooper and Boyd [1983] asymptotic analysis 
that short wave (long wave number) instability could be stabilized by surface tension. 

Thus we plot c,. vs. k as in Fig 3.6 with non dimensional surface tension 0.01 keeping all 

other parameters same. Then we find that for sufficiently high value of T and given 
value of surface tension, the high wave number modes are stabilized but some of the 
intermediate wave numbers are still unstable and very low T value does not have 
capability to stabilize the flow for all wave number. We further increased the value of 
surface tension to 0.05 and plotted in Fig. 3.8. From these Figs it is clear that for 
complete stability, we need to choose combinations of intermediate value of T and 
sufficiently high value of surface tension as presented in Fig 3.6 and 3.8. But the flow is 
still not completely stable for very low as well as very high value of T (very high or 
very low value of elasticity constant of deformable wall). Thus this needs further 
explanation why even a very high value of T is not able to stabilize flow as compared to 
intermediate values where flow is stable. For this purpose it is useful to plot neutral 
stability curve for T vs. k. Thus we plotted T vs. k in Fig. 3.17-3.24 by which it is clear 
that there are two modes of instability for the two layers Couette flow past a deformable 
wall. First one is associated with fluid-fluid interface which we call Mode-1 and other 
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and other one is associated with fluid-deformable wall interface which we named as 
Mode-2. This will be discussed in 3.4.4. 




0.001 0.01 0.1 h 1 10 100 

(a) For low value off ; c. vs. k (above) (b) For high value off : vs. k (below) 
Fig. 3.8 Variation ofc,. with k for different value off , P = 0.4, interfacial tension 
= 0.0, = 0.5, H = 0.5 and = 1 . 


54 




Fig. 3.9 Variation ofc^with k for different value ofF, P = 0.4, interfacial tension 
= 0.0, /i, = 0.5, /f = 0.5 and Re=\. 
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Fig. 3.13 Variation ofc^with k for di/lerent value ofP , p = 0.4, interfacial tension 
= 0.05, //^ = 0.5,//=0.5and7?e= 1. 

3.4.4 Fluid flow past a rigid wall is stable and deformable wall has 
destabilizing effect (//,.> 1 & P < 1/2 or /i,, < 1 <& p > 1/2) 

In this section we will discuss those cases where two-layer Couette flow past a rigid wall 
is stable (known as the thin layer effect) and deformable wall causes destabilization in the 
low k limit consistent to the equation (c). As shown in Fig. 3.14(a) flow is stable for r = 
0 expect very high wave number where it becomes unstable but from Fig. 3.14(b) it is 
clear that if we put non zero value of surface tension it becomes stable for all values of k. 
But further from Fig. 3.15 reveal that low wave number perturbations are unstable when 
r value is non-zero and extent of instability further increases as F value increases. But 
for higher wave number perturbation F has no effect on the stability of flow as the 
perturbations are confined near the fluid-fluid interface. Again in Fig. 3.15, we plotted 
vs. k which reveals that presence of surface tension would stabilize high wave number 
perturbation but small wave number perturbation can not any effect of surface tension 
and it is unstable for non-zero value of F . 
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(a) Cj vs. k: with interfacial tension = 0.0 and different values ofF . 
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(b) c, vs. k: with interfacial tension = 0.0 and different lower values ofT . 


Fig. 3. 1 5 Variation of vs. k for different values ofF , P = 0.8, interfacial tension - 0.0, 


//,, = 0.5, //=2.0 and Re=]. 
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3.45 Neutral stability curves (Tvs. k) 

In this section, we intended to find the range of f and surface tension where our 
system of flow is stable. From fig 3.17, which is a neural stability curve for parameter P = 
0.4, interfacial tension = 0.00, /u^= 0.5, H = 0.5 and Re = 1, it is clear that flow is 

unstable for any value of f . It is expected since without surface tension we can not 
stabilize short wave as discussed above. Further in fig 3.17-3.24 we have two mode of 
instability first mode is associated with fluid-fluid interface instability, which requires a 
minimum value of T to stabilize this mode and if we choose more than this value then 
mode-1 is always stable. That also reveals the reason for instability in fig. 3.9 when we 
choose gamma 0.001 which is lower than the require value to stabilize mode-1. Further 
the upper curve corresponds to mode-2 i.e. fluid-deformable wall instability, which is 
stable up to a maximum value of f beyond that flow become unstable. That’s why when 


Neutral stability curve for surface tension =0.0 



Fig. 3.17 Neutral stability curves (f vs. k) for fluid-fluid interfacial mode 
(Model) and Fluid- deformable wall mode (Mode 2) for parameter p = 0.4, interfacial 
tension = 0.00, = 0.5, H - 0.5 and Re = \. 
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Fig. 3. 1 8 Neutral stability curves (F vs. k) for Model and Mode 2 for parameter 
P = 0.4, interfacial tension = 0.01, //,. = 0.5, H= 0.5 and 7?e = 1 . 

Neutral stability curve for surface tension = 0.05 
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Fig. 3.19 Neutral stability curves (F vs. k) for Model and Mode 2 for parameter 
p = 0.4, interfacial tension = 0.05, /i:^ = 0.5, /■/= 0.5 and Re = ]. 
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Neutral stability curve for non zero surface tension 



Fig. 3.20 Neutral stability curves (F vs. k) for Model and Mode 2 for parameter 


(3 = 0.4, interfacial tension = 0.01, //,. = 0.5,/f=4 and Re - \ . 
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Fig. 3.2 1 Neutral stability curves (F vs. k) for Model and Mode 2 for parameter 
P - 0.6, interfacial tension = 0.01, = 2.0, H = 0.\ and Re= \ . 




K 

Fig. 3.22 Neutral stability curves (F vs. k) for Model and Mode 2 for parameter 
p = 0.2, interfacial tension = 0.01, ju^ = 0.5, H= 0.5 and Re = 1. 



K 

Fig. 3.23 Neutral stability curves (F vs. k) for Model and Mode 2 for parameter 
p = 0.8, interfacial tension = 0.0 1 , //,. = 2.0, /f = 0. 1 and 7?e = 1 . 
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Fig. 3.24 Neutral stability curves (F vs. k) for Model and Mode 2 for parameter 
P = 0.8, interfacial tension = 0.01, ju^ = 0.25,H=2.0 and Re= 

we plotted c. vs. k then stability of fluid firstly increases with F but after a maximum F 

limit flow becomes unstable. It is clear that high value of F certainly has ability to 
stabilize fluid-fluid interface instability but mean time it destabilizes fluid-deformable 
wall instability. Thus we need to choose an intermediate value of F and surface tension 
for which both the modes become stable. 

Above explanation is only valid for the case where flow is unstable in a rigid 
channel and deformable solid wall has stabilizing effect. Now we like to explain the case 
when flow is stable in rigid channel and we intended to destabilize with the help of 
deformable solid wall. From Fig 3.24 it is clear that low wave number perturbations at 
the fluid-fluid interface are unstable if we increase the F value beyond the maximum 
limit and Fig 3.14 and 3.15 reveals that surface tension has stabilizing effect on the short 
waves. 
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3.4.6 Effect of various parameters on the stability 

(a) Effect of viscosity ratio on the stability of flow system; 

First of all we want to know the effect of viscosity ratio ) on the stability of two-layer 

fluid flow past a deformable wall for that purpose we plotted F vs. k in Fig. 3.25, which 
reveals that as the value of increases system needs lesser value of F for flow to 
become unstable. We can explain this observation by considering that the non- 
dimensional velocity gradient at the solid-deformable wall interface, which is given by 
jiij{\ + l3 (//^ -■!)), and for a given value of F, this quantity increases with increase in 
ratio of viscosities, thus having a destabilizing effect on the fluid-solid interfacial mode. 
Similarly, for decreasing ratio of viscosities, the non-dimensional velocity gradient at the 
fluid-solid interface decreases, thus having a stabilizing effect. 



Fig. 3.25 Neutral stability curves (F vs. k) for Mode 2 for the parameter P = 0.4, 
interfacial tension = 0.0, //= 5.0 and Re= \ . 

(b). Effect of density of upper fluid on the stability of system of flow: 

Although density difference itself has instability modes (density stratification modes) but 
here we intended to examine the effect of density difference on stability of fluid-fluid 
interface. From Fig. 3.28 it is clear that as density of upper tluid increases as compare to 
density of lower fluid, flow become more and more unstable. 
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Fig. 3.26 Variation of F vs. k for different values of densities ratio, |3 = 0.6, 
interfacial tension = 0.01, //,.= 0.5, H=2.0 and Re=l. 

(c) Effect of height of deformable wall on the stability of system of flow: 



Fig. 3.27 Neutral stability curves (F vs. k) for Mode 2 for the parameter |3 = 0.4, 


interfacial tension = 0.0, H= 5.0 and Re=l. 
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3.5 Conclusion and future work 

We have shown in this chapter that the deformable wall plays an important role in 
the stability of two-layer fluid flows. This can be explained on the basis of coupling 
between two interfacial perturbations. But it has certain limitations, as it is only able to 
affect long waves and finite waves, since coupling between two interfacial perturbations 
is not possible for short waves. This is because, the velocity perturbations are confined 
near the interfaces only. Further short waves are stabilized by the interfacial tension 
between the two fluids. 

Here we intended to analyze the instability due to viscosity stratification. For this 
purpose we first considered the case where the interface is unstable in a rigid channel and 
deformable wall has stabilizing effect. But we had shown that it is not possible to 
stabilize the two-layer flow by using any deformable wall rather it depends of the 
coefficient of elasticity of deformable wall. We found a region of F for which our 
system of flow becomes stable. But if we have a system such that the F value is less 
than of minimum transition value F then fluid-fluid interface becomes , unstable and 
simultaneously when we have F value more than maximum transition value then fluid- 
deformable wall interface becomes unstable. Next we consider a case where two-layer 
flows are stable (thin layer effect) and deformable wall has destabilizing effect. Similar to 
the above case in this case also there is a region of F which deformable wall is able to 
destabilize the system of flow. Finally, we had shown the effect of various parameters on 
the stability of fluid flow. For this, first we choose a case in which fluid-fluid interface is 
stable and (luid-deformable wall stability increases as increases. Then we had shown 

that as Reynolds increases and height of deformable wall have destabilizing and 
stabilizing effect as their values increases. Finally we had shown that the density 
difference also has effect on the both fluid-fluid mode as well as fluid deformable wall 
modes. 

Although we talk about the coupling between two interfacial perturbations, we 
never got any information that how actually they are coupled with each other for that 
purpose we mu.st need to know the perturbation growth for this system of flow. Other 
than that in present analysis we predicted stability and instability based on linear stability 


70 



analysis but we can further conforms the results by doing non-linear stability analysis. 
We also showed that density difference have effect on instability caused by viscosity 
stratifications. There is also the effect of viscosity difference on the in the stability 
caused by density stratification. In fact these things are already available in literature for 
rigid wall case but now it is possible to study these effects for the case of deformable 
wall. This present study is associated with finite thick fluid layers so if we reduce the 
thickness of (luid layers to the range of thin layer where vander waals force are 
significant. 1'hen our analysis is also valid for the case of thin film but we have to 
incorporate tlie vander waals forces. 
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CHAPTER 4 

Conclusion and future work 

We found that Couette flow past a deformable wall has an inviscid mode, which is very 
important for these types of flows. If one tries to study flow past a deformable wall 
experimentally then one will encounter inviscid modes, because they require lower value 
of r . Further we confirmed that Poiseuille flow also have both wall modes, as well as 
inviscid modes. However there is a difference between these two cases. In Couette flow 
one of the wall modes captures the inviscid mode when we extend wall mode solution 
numerically. But in I'oiseuille flow, instead of wall modes, zero Reynolds number modes 
capture the iiwiscid mode. In case of Poiseuille flow, we got the same scaling as Couette 
flow but velocity field distributions have different nature instead of expected critical layer 
for inviscid modes. The present study needs further investigation of velocity distribution 
in the case of I’oiseuille flow. Other than this, our study numerically predicted inviscid 
modes (or both the plane Couette as well as plane Poiseuille flow. The prediction of 
inviscid modes for plane Couette flow is in stark contrast with the asymptotic prediction 
of Kumaran 1 1995], who concluded that the inviscid modes are stable in the case of plane 
Couette flow. Hence we further need asymptotic explanation for the inviscid cases. 

We have shown in Chapter 3 that the deformable wall plays an important role in 
the stability of two-layer fluid flows. This can be explained on the basis of coupling 
between two interlacial perturbations. First we have shown that deformable wall is able 
to stabilize those flows, which are unstable in rigid channels. In the other case where flow 
is stable in rigid channel, deformable wall has destabilizing effect on the two-layer 
interlacial modes. Finally, even different density fluid flow can be stabilized by the 
deformable wall. 

Although we talk about the coupling between two interfacial perturbations, we 
never obtained any information that how actually they are coupled with each other for 
that purpose we must need to know the perturbation growth for this system of flow. Other 
than that in present analysis we predicted stability and instability based on linear stability 
analysis but we can further confirm the results by doing non-linear stability analysis. We 
also showed that density difference has an effect on instability caused by viscosity 
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stratifications. There is also the effect of viscosity difference on the in the stability 
caused by density stratification. In fact these things are already available in literature for 
rigid wall case but now it is possible to study these effects for the case of deformable 
wall. This present study is associated with finite thick fluid layers so if we reduce the 
thickness of fluid layers to the range of thin layer where van der Waals force are 
significant. Then our analysis is also valid for the case of thin film but we have to 
incorporate the van der Waals forces. 
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